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Abstract 

The water uptake by roots of plants is examined for an ideal situation, with an 
approximation that resembles plants growing in pots, meaning that the total 
soil volume is fixed. We propose a coupled water uptake-root growth model. A 
one-dimensional model for water flux and water uptake by a root system growing 
uniformly distributed in the soil is presented, and the Van Genuchten model for 
the transport of water in soil is used. The governing equations are represented 
by a moving boundary model for which the root length, as a function of time, 
is prescribed. The solution of the model is obtained by front-fixing and finite 
element methods. Model predictions for water uptake by a same plant growing 
in loam, silt and clay soils are obtained and compared. A sensitivity analysis to 
determine relative effects on water uptake when system parameters are changed 
is also presented and shows that the model and numerical method proposed 
are more sensitive to the root growth rate than to the rest of the parameters. 
This sensitivity decreases along time, reaching its maximum at thirty days. A 
comparison of this model with a fixed boundary model with and without root 
growth is also made. The results show qualitative differences from the beginning 
of the simulations, and quantitative differences after ten days of simulations. 
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1. Introduction 


In the development of a theory to describe plant water uptake, electrical 
analogues of the system have been used for analysis The analogues 

are based on the assumption that rooting patterns are uniform and constant in 
each soil layer. Steady flow is presumed in both the soil and the plant over the 
period of calculation. In this approach the plant water potentials are primarily 
the result of an imposed value of transpiration rate and its variations. Later 
papers iili have presented detailed reviews on plant water uptake. In those 
papers the Richards equation is used, with a sink term. Another approach is to 
model the water movement and uptake over large areas, using individual plant 
0 , or global behavior Q . A microscopical approach has also been proposed Q , 
where the total water uptake is calculated based on using a constant value for 

m\ — 


the entire rooting profile. In more recent papers [lOj, I 111 . Il2| the root growth has 


been taken into account, still using a fixed domain. The root growth is inscribed 
on a domain that is not a function of time. Some other papers about nutrient 
uptake consider root growth and instantaneous coupling with the nutrient flux 
by using a variable domain approximation 131. In this last model [l^ a variable 
root length, and consequently, a variable available volume of soil to each root of 
a root system is considered using a moving boundary model. In this model the 
root system is uniformly distributed in the soil and the variation of available soil 
volume per unit of root length is modeled by a moving boundary. The approach 
presented here is based on that in [l^ . In the proposed model plants growing 
in controlled conditions, as in a growth chamber, are assumed. A constant 
temperature and evapotranspiration rate is presumed. In this situation, the 
water potential at the root surface is determined by the soil water potential, 
and consequently determines water uptake by the growing root system. The 
proposed model considers an uniform root water uptake for all the root system. 
The goal of this paper is to present a simplified model of water uptake coupled 
with a growing root system and analyze the influence of system parameters on 
water uptake using typical values. 


2. Model 


Darcy’s law describes the flow of water on a porous unsaturated medium as 

J(r, t) = — A'('I'(r, t))V'I'(r, t) (1) 

with J[cm^/crn^s] the water flux per surface unit at position f\cm] at time t[s], 
K[cm/s] the soil water conductivity and 4'[cm] the soil water potential. 

The corresponding continuity equation (mass conservation) is given by 


-V-J = 




( 2 ) 


2 







t = 0 


t = t 


Figure 1; Homogeneous rooting in soil and its time evolution. 


with 9[crrt’/cm^] the soil water content per unit of volume, with the approxi¬ 
mation of only radial flux the transport equation results 


5 rx / M 1 ^ 


rK{^{r,t))— [^'(r,f)] 


( 3 ) 


where r is the cylindrical radial coordinate, and 


C'(^'(r,t)) = 


( 4 ) 


is the differential capacity of water [cm~^]. 

The soil water constitutive relations for K, 6 , and C, as functions of 4', are 
the ones proposed by Van Genuchten [1^, and consist of the expressions given 
by: 



where Ks[cm/s] is the saturated soil conductivity, 9s[cm^/cm^] is the saturated 
soil water content, is the residual water content, 'Pe[cTO], p[l] and 

n[l] are experimental coefficients, and m = 1 — 1/n. 

We presume that the water potential does not change with soil depth, the 
soil does not evaporate, and laboratory conditions, like light and temperature 
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are maintained constant. The root density is homogeneous on the soil, the total 
volume is fixed (as in pots), therefore the soil volume per unit of root length is 
decreasing as in Figure [TJ Based on these assumptions, and taking into account 
the root length density as a function of time (t) the following moving boundary 
model in cylindrical coordinates @,[11 is proposed 


«'(r,0) 

G{^{so,t)) 

m 


1 d 
rC(4'(r, t)) dr 
(j){r) 




d'i' 

-2TTSoK{'^{so,t)) — {so,t) 



( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 


with So < r < R{t) and 0 < t < T^ax, where T^ax is the maximum time for 
which the system has meaning {R{Tmax) > sq), sq is the root radius, Rq is the 
initial half-distance among roots, R{t) is the instantaneous half-distance among 
roots (a decreasing function as root density grows), Iq is the initial root length, 
and l{t) is the instantaneous root length. Equation ([8]) is the pressure head based 
Buckingham-Richards equation. The condition (jOj is the initial water potential 
profile, with (j){r) a single valued function. The condition (|10ll represents the 
flux (H) on the moving boundary r = R{t), which will be considered null in this 
paper as an approximation to a soil isolated. The condition (EH) is the boundary 
condition at the root soil interface (^ = 50 ) representing the root water uptake 
per unit of root length (^G{'i’{so,t))[cm^/cm s]). For the water uptake function 
(G'('I')) we use the function proposed by Feddes d which is given by: 


Gi'f) =0 if 

G^^-) = Smax if 

Gi'f) = S^axi^^ if 
G(^') = 0 if 


0 > vk > 

vki > vk > VI/2 
'k2 > 'h > 'f'3 

^3 > ^ 


where 'I'l, '^2 and ^'3 are the anaerobiosis point, the limiting point, and the 
wilting point, respectively. Smax [cm^/s] is the maximum water uptake per 
unit of root length. A diagram of this function can be seen on Figure [2j The 
condition EH is the time dependence of the moving boundary {R{t)) which is 
obtained presuming a fixed total volume including soil and roots, and a linear 
growth rate {l{t) = Iq + V t) [l^, where V[cm/s] is the root length growth 
rate. Graphical evolution of this system with time can be seen in Figure [TJ A 
schematic mathematical diagram of the problem is shown in Figure |3| Once 
the water potential ^'(so,^) at the root surface is obtained, the water uptake is 
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Figure 2: Root water uptake function as proposed by Feddes [J|. The zone (a) is the zone of 
maximum water uptake, while the zone (b) is the zone of linear water uptake. 


computed with a variable domain integration method given by [if 


U{t) = lo f G{'i'{so,T))dT 
Jo 


+ 


G{'il{so,T))dT i{T)dT 


(13) 


where U [cm^] is the cumulative water uptake at time t. This last expression 
can be simplified to (see [Appendix AD 


U{t) = [ G{^{so,T))liT)dT, 

Jo 


then the instantaneous water uptake can be defined as 
t7,(t) = t7(t) = G(^(so,t))Z(t). 


(14) 


(15) 


The remaining soil water {W [cm^]) as a function of time can be calculated 
using 


rm 

W{t) = 2Trl{t) / 0 (4'(r, t)) r dr, 

J sn 


(16) 


if the root is not growing then l{t) = lo = constant, and as a consequence 
R(t) = Rq, therefore the water remaining in the soil in this case is 


rRo 

W{t)=2TTlQ / 0 (4'(r, f)) r dr. 

J Sn 


(17) 
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Figure 3: model of the domain of validity of the proposed model, the shaded zone (a) represents 
the zone where the equation ((Sj is valid. The line (b) is the initial value condition |[9]l. In the 
curve (c) r = R{t) the boundary condition IIIQII is used. In the line (d) the root water uptake 
condition {HI is used. 


To solve the system ®-(ini) the domain is transformed to a dimensionless 
form by using the following expressions: 


X = 


^{x,t) = 

or their inverses given by 


r- So 
R{t) - So ’ 
t 

to ’ 

5'(x(i?(t) - So) + so,Tto) 






r = So + {R{t) - so)x 
t = toT 

'h(r, t) = (a;,T) ='he'h 


r - So ^ 

R{t) - So' to) ' 


(18) 


(19) 


6 









and the following definitions 


k{^{x,T)) = 


(20) 

Ks [1 + $"(x,t)]"^(^+") 

n 

II 


(21) 

{9s-9R){l-n) [1 + 

v{t) = 

R{t) 

50 

(22) 

w{t) = 

m 

So 

(23) 

II 

G(vE-(r,t)) 

(24) 

H{'it{r,t)) = 

0, 

(25) 

c*{<^{x,t),x,t) = 

((u(t) - l)x + 1) c(<i)(x, r)) 

(26) 

k*{^{ x,t), x,t) = 

((u(t) - l)x + 1) k{<^{x, t)) 

(27) 

c($(a:, T),x,r) = 

xv{t)w'{t) ^ 

— 7 :—, w / X- ttC ($(a;,r),a:,r) 

2n;(r)(u(T) - 1) ^ ^ ^ 

(28) 

v{t) = 

Ro 1 ^0 

(29) 

So V sow{t) 

7 = 

to'l’eKs 

(30) 

{9s - 9ii){l - n)sl 

a = 

Smax 

(31) 

2TTKs'^e 

'^max — 

T 

max 

io 

(32) 


where k{^{x,T)) and c(<i)(a;,r)) are the dimensionless forms of and 

t)) respectively. The dimensionless number 7 takes into account soil 
properties and the geometrical proportions of the root, to is chosen to make 
7 = 1. The dimensionless number a takes into account one flux per unit of 
length, related to the water uptake, the term Ks'^e is a flux per unit of length 
when the soil is saturated and the gradient of the matric potential is equal to 
dte/lcm. 

Taking into account (fTSl) and the definitions (l20l) - (l32]) . the system (l8l)" (fT^ 
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is transformed in a dimensionless form in the domain (0,1) x (0, Tmax) given by: 


c*($(a;,r),a;,T)— [$(a:,r)] 


$(a:, 0) 




cr g($(0,T)) 


c{^{x,t),x,t)— [$(a;,T)] - 


A 

dx 
ip{x), 


7 fc*($(a;, r), cr, r) d 
{v{t) — 1)2 dx 




(33) 


(34) 


0 , 


(35) 


.;;i^e(4.(0,r),0,r)g(0,T), 


(36) 


where i^{x) = (^(sq + (-Ro ~ so)2^)/^e, is the dimensionless initial profile. This 
transformation maps the spatial domain variable in time to a fixed domain in 
time, and adds a term on the right side of the transformed transport equation 
(1331) which contains the variation of the moving boundary. In order to solve 
the model (I55H5B1) the non-linear finite element method [id is applied and the 
resulting model is solved by using the software FlexPDE [l7| with an adaptive 
mesh of around 400 nodes. 


3. Results 


All simulations were performed by a same hypothetical plant in three types 
of soils (loam, silt and clay) for the same soil-root volume (4000cm^) with the 
same total water content {6i = 0.30). This initial condition is fixed using the 
’’field capacity” concept given by Ritchie [l^. Hydraulic soil data selected were 
those for loam, silt, and clay based on Q. The soil parameters used are shown 
in Table m For the plant parameters different sources were used. The values 
of 4 '2 and 4'3 were taken from [l^, the value of Ti was chosen to assure that 
wateruptake was possible at the beginning of the simulation, sq and V are 
typical values from the literature , and Iq was chosen to be 1 cm to simulate 
a plant at the start of growth. The plant and soil volume parameters are listed 
in Table [2] Table |3] shows parameters of soil and plant properties and the initial 
“available water” (t/g)- This last parameter is approximated as 


Uo = [di - 0(4-3)] X 4000cm^ 


(37) 


and represents how much water can be extracted from the soil before it reaches 
an uniform water potential at the wilting point. Beyond this point the root 
cannot extract more water. 

Figures m [5l and [6| show the soil water potential profiles at different times, 
with an initial moisture condition equivalent (0 = 0.30), for the loam, silt and 
clay soil, respectively. The curves reveal that a high water potential gradient is 
developed in a very small time period for the clay soil, while for the loam and silt 
soil the development of the water potential gradient is more gradual. For loam 
and silt soil the root dries the soil near the root surface on the first 30 days. After 
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Soil 

Os 

0R 

Ks (cm/s) 

5'e (cm) 

n 

P 

Loam 

0.43 

0.078 

2.89 X 10-4 

-27.78 

1.56 

0.5 

Silt 

0.46 

0.034 

7.00 X 10-5 

-62.5 

1.37 

0.5 

Clay 

0.38 

0.068 

5.80 X 10"® 

-125 

1.09 

0.5 


Table 1: Soil hydraulic properties 


Parameter 

Value 

^0 

1 

cm 

i?o 

35.7 

cm 

V 

1 * 10-* 

cm/s 

So 

0.05 

cm 

Synax 

2 *10-® 

cm^/s 

4-1 

-1 

cm 

4'2 

-750 

cm 

4'3 

-17500 

cm 


Table 2: Values for the plant and soil volume parameters 


Soil 

'I'(r, 0) (cm) 

to (s) 

Uo (cm®) 

Loam 

-51.4 

0.0614 

850 

Silt 

-194 

0.0901 

852 

Clay 

-3280 

0.0968 

128 


Table 3: Simulation parameters depending on soil and plant properties 
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Figure 4: Water potential profiles at the scaled domain for the loam soil. 



Figure 5: Water potential profiles at the scaled domain for the silt soil. 


that period these soils show a water potential on the root surface close to the 
wilting point. The root starts to retrieve water near the moving boundary R{t). 
This is shown by an increase of the modulus of the water potential (decrease of 
his value) at x = 1 for the dimensionless domain, or r = R{t) for the physical 
domain. The clay soil shows a similar behavior but the time period at which 
the potential near root zone reaches potentials close to the wilting point is 10 
days. 

Figure [7] shows the relative water uptake (G'('I'(so, t))/Smax) per unit of root 
length as a function of time for the three types of soil, at different root growth 
rate. At the initial time the water uptake for loam and silt soil develops at the 
zone of maximum uptake (zone (a) in Figure [2]), while clay soil does the same 
for the linear uptake zone (zone (b) in Figure [2|). The figure shows that the 
water uptake does not vary much for the first 20 days as a function of V for 


10 





































Figure 6: Water potential profiles at the scaled domain for the clay soil. 


loam soil, after that period loam soil enters on the linear regime at different 
times depending on the value of V. After 40 days the gap among the curves 
on the log scale remains almost constant and all the curves are in the linear 
uptake regime. This means that the instantaneous root water uptake curves 
have a multiplicative constant among them. Silt soil has a very similar behavior 
to loam soil but the linear regime occurs at shorter times. For clay soil the 
period of similar water uptake is 10 days and after that period the curves begin 
a gradual separation that continues until the end of the simulation. Clay soil 
also shows a very sharp initial decrease on the curve (one order of magnitude), 
this is due to the linear regime and a large decrease of 4'(so,i) at the beginning 
of the simulation. 

To compare the above simulations with those for a constant root length 
density in a fixed domain model a simulation with V = 0 {R(t) = Rq) was made. 
For the comparison to be useful the root length was taken to be the average root 
length {Iq = I = 589.8cm). Using this value the domain is fixed and the moving 
boundary model then becomes the model proposed by Personne @ . The influxes 
on the root surface estimated by the fixed boundary model are then integrated 
using d with the fixed root length, to compute the cumulative uptake. This 
model will be referred to as Fixed Boundary Fixed Length (FBFL). Similar to 
published results on nutrient uptake (e.g. Claasen and Barber [21], Cushman 
|22|) the influxes obtained with the fixed boundary model can be integrated 
using equation m with a variable root length, to compute the cumulative 
uptake. This model will be referred to as Fixed Boundary Variable Length 
(FBVL). Figure m shows the instantaneous root system water uptake against 
time (i.e. G(4'(so,t)) x l{t)), for the three models. The time at which the 
instantaneous water uptake is maximum is called the Maximum Uptake time 
(MUT). The MUT is only present when the roots are growing. The straight line 
in the beginning is caused by the assumption of linear root length growth and 
the constant water uptake (G('L) = Smax), on the models with root growth. In 
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Figure 7: Relative water uptake per unit of root length time evolution for the different soils. 
The curves refer to changes in root growth speed (V). The square is a zoom on the first day 
of simulation for the clay soil. 




contrast when the root is taking up water in the water stressed (linear zone) of 
the water uptake function the instantaneous water uptake as a function of time 
is non linear, in all models. We also observe that in the clay soil for the models 
with root growth the water uptake initially decreases given the large decrease 
in soil water potential at the root surface. This variation is shown in the inset, 
the later increase in water uptake is due to the root growth. The models with 
growing roots have a similar course in time, but the differences between the 
MUT-values are consistently showing that the MUT in the FBVL occurs earlier 
after the start of root growth. The instantaneous water uptake of the FBFL 
model and the FBVL models differs only by the multiplicative factor 

Figure IHl shows the cumulative water uptake as a function of time for the 
three models. The cumulative water uptake at 90 days (17(90)) is close to the 
amount of water initially available {Uq) for the FBFL and the moving boundary 
models (see Table O. For the FBVL model the cumulative water uptake is lower 
than that predicted by the other. This difference is due to the smaller value 
of the MUT. After 90 days 20% of the initially available water in the clay soil 
is remaining, whereas for the silt and loam soil it is in the order of 10%. The 
clay soil does not trend to a constant value due to the less water uptake. It 
is expected that over a longer period the moisture content in the clay soil will 
asymptotically approach to a lower value, that might be on the proximity of the 
values found on the other soils. 

To analyze the effect of parameter variation on model output variables, sev¬ 
eral sensitivity diagrams were made using a local approach [23|. The results 


are presented in Tables and [3 as a function of the multiplicative fac¬ 

tor of each parameter for the loam soil. The more relevant parameters are V 
and Smax, having both positive correlation on the water uptake at 30, 60 and 
90 days and negative correlation with the MUT. As for the WUMUT, Smax 
has a negative correlation and V has a positive correlation, sq has an almost 
constant sensitivity and always positive correlation. Therefore the experimen¬ 
tal measurement errors on those parameters would be more amplified on the 
output of the model. For the silt soil there is a very similar pattern than for 
the loam soil. For the clay soil there are some differences. Smax is no longer a 
relevant parameter and sq is a relevant parameter with positive correlation to 
the accumulative water uptake. For the MUT and WUMUT sq has a negative 
correlation. The behavior of V is left unchanged with respect to the loam soil. 
MUT shows almost no sensitivity to Smax but the WUMUT shows that this 
parameter is very relevant with positive correlation. All the sensitivities are non 
linear with some exceptions for sq- 

Mass conservation implies that the total water volume must be the same at 
the beginning and at an arbitrary time of the simulation. Since the water can 
only be in the soil or inside the root, therefore the total water volume (T [cm^]) 
at time t would be 

T{t) = U{t) + W{t) (38) 

with U{t) and W{t) defined as in equations (fldll and (IT6l) respectively. At the 
beginning of the simulation the cumulative water uptake is null, therefore the 
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Fixed Boundary Fixed Length (FBFL) 



Fixed Boundary Variabie Length (FBVL) 



Moving Boundary 



Time [days] 


Figure 8: Instantaneous water uptake time evolution for the three models to be compared, 
from top to bottom fixed boundary without root growth, fixed boundary with root growth 
and moving boundary by the root growth. 
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Fixed Boundary Fixed Length (FBFL) 



Fixed Boundary Variable Length (FBVL) 




0 10 20 30 40 50 60 70 80 90 


Time [days] 


Figure 9: Cumulative water uptake time evolution for the three models, respectively fixed 
boundary without root growth, fixed boundary with root growth and moving boundary with 
the root growth. 
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Factor 

0.5 

1.5 

2.0 

U30 

Ueo 

U90 

MUT 

WUMUT 

-1.5 10"^ 
-1.6 10-4 
-9.1 10-5 
2.1 10-5 
-7.8 10-5 

1.5 10-5 

1.6 10-4 

8.7 10-5 
-2.1 10-5 

4.6 10-5 

2.9 10-5 
3.2 10-4 
1.7 10-4 
-4.2 10-5 
1.1 10-4 


Table 4: Sensitivity of the moving boundary model to the variation of the initial length Iq 
for the loam soil. C/30, C/eo, s,nd C/90 ^-^e the cumulative water uptake at 30, 60 and 90 days 
respectively. MUT and WUMUT are the maximum uptake time and the cumulative water 
uptake at maximum uptake time respectively. The sensitivity is the change of the output 
divided by the value of the output without variation of the parameter. 


Factor 

0.5 

1.5 

2.0 

U30 

-1.5 10-^ 

9.7 10-5 

1.7 10-2 

Uqo 

-1.0 10-2 

6.6 10-5 

1.2 10-2 

Ugo 

-7.7 10-5 

4.7 10-5 

8.1 10-5 

MUT 

-1.4 10-2 

9.0 10-5 

1.6 10-2 

WUMUT 

-2.8 10-2 

1.8 10-2 

3.2 10-2 


Table 5: Sensitivity of the moving boundary model to the variation of the root radius 5o for 
the loam soil. C/30, Ugo are the cumulative water uptake at 30, 60 and 90 days 

respectively. MUT and WUMUT are the maximum uptake time and the cumulative water 
uptake at maximum uptake time respectively. The sensitivity is the change of the output 
divided by the value of the output without variation of the parameter. 


Factor 

0.5 

1.5 

2.0 

U30 

-4.5 10-5 

6.4 10-2 

8.0 10-2 

Ueo 

-1.3 10-2 

1.3 10-5 

1.5 10-5 

U90 

-1.1 10-5 

1.0 10-5 

2.3 10-5 

MUT 

4.8 10-4 

-2.1 10-4 

-3.4 10-4 

WUMUT 

8.7 10-2 

-6.3 10-2 

-1.1 10-4 


Table 6: Sensitivity of the moving boundary model to the variation of the maximum root water 
uptake per unit of root length Smax for the loam soil. C/30, C/eo, and C/90 are the cumulative 
water uptake at 30, 60 and 90 days respectively. MUT and WUMUT are the maximum 
uptake time and the cumulative water uptake at maximum uptake time respectively. The 
sensitivity is the change of the output divided by the value of the output without variation of 
the parameter. 
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Factor 

0.5 

1.5 

2.0 

U 30 

-4.5 10"^ 

1.0 10-1 

1.5 10-1 

Uqo 

-5.9 10-2 

2.5 10-2 

3.9 10-2 

U 90 

-3.6 10-2 

1.6 10-2 

2.6 10-2 

MUT 

4.1 10-1 

-1.8 10-1 

-2.9 10-1 

WUMUT 

-7.3 10-3 

4.3 10-3 

7.4 10-3 


Table 7: Sensitivity of the moving boundary model to the variation of the maximum root water 
uptake per unit of root length Smax for the loam soil. U 30 , Ueo, and U 90 are the cumulative 
water uptake at 30, 60 and 90 days respectively. MUT and WUMUT are the maximum 
uptake time and the cumulative water uptake at maximum uptake time respectively. The 
sensitivity is the change of the output divided by the value of the output without variation of 
the parameter. 


total water is the water in soil. With the above considerations the change of the 
cumulative water uptake plus the water remaining in the soil minus the initial 
water content {W{t) + U{t) — W{0)) as a function of time was calculated. For the 
FBVL model the mass balance cannot be done because to compute the water 
remaining in soil, the operation must be done with pressure head profiles as a 
function of time which has been calculated in fixed domain, but this result must 
be compared with the cumulative uptake by a growing root, which has been 
calculated integrating in a variable domain. The results of those calculations 
are shown in Figure [TOl It shows that there is an effect of water mass loss for the 
moving boundary model that is not present on the fixed boundary model. There 
are two possible causes for this mass loss. The first cause is the assumption of 
total volume constant in the formulation of the model, reflected in the moving 

boundary formula (^R(t) = i.e., the volume that is kept constant is 

the soil plus root volume, being the volume occupied by the root is Vr = Trs^lit). 
Then the soil volume is the total volume minus the root volume, therefore the 
root occupies a bigger fraction of the total volnme as time increases (see Figure 
[J). The water content in the volume of soil removed by the root is not taken 
into acconnt in the calculation of the water in the soil volnme {W{t)). This 
means, in other words, that, for the model, the root is “eating” the soil, and 
the water which is in that portion of the soil. The second cause is the numerical 
errors. This case is similar for both models and does not contribute a large mass 
loss. An important remark is that the FEM is used to solve the water potential 
('k), not the water content (0), on a discretized space with a finite precision, 
and cumulative errors could lead to a mass creation or loss depending of the 
parameters. Figure (TU] shows the soil volume loss effect plus the water mass loss 
due to the used numerical method in each integration step. In each soil for the 
moving boundary model the water loss is around 1% of the total water. For the 
clay soil the total loss is around 15% of the total water uptake, while for the 
other soils the loss is less than 2%. The water loss for the fixed boundary model 
is negligible compared with the total water or the cumulative water uptake. 
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Fixed Boundary Fixed Length (FBFL) 




Figure 10: Water mass balance {W{t) — W(0)) versus time (t) for the values shown in 

Table [2 for the FBFL and the Moving Boundary models. The total water is the cumulative 
water uptake plus the water remaining in the soil, the initial water is the total water at the 
start of the simulation. 
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4. Conclusions 


An important remark is that all the results presented before are computed 
for a single plant, in a fixed soil volume with initial soil water content equal for 
all simulations, to compare the effect of the different soils on the water uptake. 
The parameters used represent traits of the plant. V can be regarded as the 
response to nutrient and water uptake and is influenced by the plant genetics. 
Smax is linked to atmospheric factors and to genetic variation, sq is the constant 
root radius. Iq is a growing stage trait, changes in it only means a change in 
the development of the rooting system when the simulation begins. 

There is a large difference in behaviour among the silt and loam soil on one 
hand and the clay soil on the other hand. Initial pressure head (Table IS]) and 
pressure head profiles are different (Figures 21 [S] and [B]), and instantaneous water 
uptake is much smaller (Figures [7] and [5]). This can be explained in terms of 
the water uptake function (G(4')). For the clay soil the water uptake function 
is on the linear regime (zone (b) of Figure [^l and the water potential is sharply 
decreasing (drying the soil) at the root-soil interface. Once the time variation 
of the water potential at the interface is stabilized the root growth has the 
dominant effect on uptake and the curve is similar to the other curves, although 
it shows a non linear behavior. 

From Figure [TOl the mass balance is not exactly zero owing to the numerical 
method used. For the FBFL model the differences are due to numerical errors 
only and are very low compared with the water volume loss on the moving 
boundary model. The mass variation is produced on the high instantaneous 
uptake zone, this is a known flaw of the finite elements method. The moving 
boundary model has a higher mass variation compared with the FBFL model by 
the accumulation of two factors, numerical errors as in FBFL model, and a soil 
volume loss effect. The soil volume loss effect is introduced by the assumption 
of constant total volume (roots grows at the expense of a decrease in the soil 
volume), which is reflected on the formula of the moving boundary (1121) . For the 
moving boundary model the differences on the total values of water mass loss by 
a soil volume loss effect and numerical errors are due to the different variations 
of the soil water contents. The clay soil has low available water (see Table El), 
but the total water mass loss is similar to the other soils, and it represents 
about 10% of the initial total water. The assumptions that generate the soil 
volume loss effect should be revised when soils for which the initial potential 
is on the linear zone (i.e., low water availability) are studied. The water mass 
loss is consistent with the model assumptions, and the numerical induced errors 
are low compared to the water uptakes (< 15% for clay soils, and < 2% for the 
other soils). There are two main actions to avoid the volume loss effect. The 
hrst one is to keep the soil volume constant and not the root-soil volume. The 
second one is to evaluate the mass loss and put it back into the soil using a 
source function or a modification of the boundary conditions. Each procedure 
will have its advantages and disadvantages, but, since the water loss is very low 
compared with the water uptake for the studied cases, the revision is left for a 
specific work on low initial available water. 
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The obtained results shows a global consistency with the structure of soils 
studied, becoming a valuable tool to study the water uptake in a more complex 
situation (for example when the effect of a variable evapotranspiration on Smax 
is considered). 

Obviously, the results would change when plants not growing in pots but in 
the held (in this case equation m is invalid) are considered. Here is necessary to 
develop a new rooting development function R(t). One change in this model to 
be useful in a held situation is the use of a coordinate representing the depth and 
water transport by gravity. Here R(t) should incorporate the root architecture. 
Moreover, since this is a hrst approach a simple water uptake function has been 
used. This function could be changed to a more complex function 2J, [25 1. 
Moreover the effect of V and Smax could change substantially after the hrst 30 
days in held conditions, since there will be circadian and climatic effects that 
will alters the root water uptake function (G), and the evapotranspiration will 
be a parameter to be measured or calculated before the simulation. 

The diherences among models with and without root growth can bee seen in 
hgures |5| and [H The dynamics of water uptake for the hxed boundary models 
are very different from those of a moving boundary model. The maximum in¬ 
stantaneous water uptake is achieved much faster for the hxed boundary model. 
Also a difference in the beginning and end of the simulation for the values of 
the instantaneous water uptake can be observed, while for the FBFL and the 
moving boundary model the instantaneous water uptake reaches values close to 
zero as time goes by, is not clear that the FBVL model has the same behavior. 
This difference is because the inhux on root surface is computed for a single root 
in a hxed domain and the cumulative uptake by the growing root is calculated 
using these inhuxes. From a physical point of view, to compute the cumula¬ 
tive uptake inhuxes calculated in a variable domain should be used. Therefore 
the mass balance cannot be calculated for the FBVL. For those considerations 
of mass conservation the FBVL approach should be avoided on water uptake 
models. 
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Appendix A. Deduction of simplified total water uptake 

The formula for water uptake given by [l5| is 
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the second integral on the right side can be reduced, by a parts integration, into 

U{t) = Iq [ G(4'(so,T))dT 
Jo 

+ (y G(^'(so,r))dT^/(r) -J (-G(^'(so,r)))/(T)dr, (A.2) 

U{t) = lo [ G{'^{so,T))dT 
Jo 

f G(d'(so,r))(ir^ ;(0) + / G(^'(so,r))Z(r)(ir (A.3) 

\Jo J Jo 


Therefore the final result is 


U{t) = f Gi^{so,T))liT)dT 
Jo 


(A.4) 
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